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Abstract Strong {B 3> 10^ G) and superstrong {B > 
10^^ G) magnetic fields profoundly affect many ther- 
modynamic and kinetic characteristics of dense plasmas 
in neutron star envelopes. In particular, they produce 
strongly anisotropic thermal conductivity in the neutron 
star crust and modify the equation of state and radia- 
tive opacities in the atmosphere, which are major in- 
gredients of the cooling theory and spectral atmosphere 
models. As a result, both the radiation spectrum and the 
thermal luminosity of a neutron star can be affected by 
the magnetic field. We briefly review these effects and 
demonstrate the influence of magnetic fleld strength on 
the thermal structure of an isolated neutron star, putting 
emphasis on the differences brought about by the super- 
strong fields and high temperatures of magnetars. For 
the latter objects, it is important to take proper account 
of a combined effect of the magnetic field on thermal 
conduction and neutrino emission at densities p > 10^" 
g cm^"^. We show that the neutrino emission puts a in- 
dependent upper limit on the effective surface tempera- 
ture of a cooling neutron star. 
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1 Introduction 

Thermal emission from neutron stars can be used to 
measure the magnetic field, temperature, and compo- 
sition of neutron-star envelopes, and to constrain the 
properties of matter under extreme conditions (see, e.g., 
lYakovlev fc Pethick|[200l lYakovlev et al.ll2005l and ref- 
erences therein). To achieve these goals, one should use 
reliable models of the atmosphere or condensed surface, 
where the thermal spectrum is formed, and of deeper 
layers, which provide thermal insulation of hot stel- 
lar interiors. In these layers, the effects of strong mag- 
netic fields can be important. In recent years, significant 
progress has been achieved in the theoretical description 
of neutron-star envelopes with strong magnetic fields, 
but new challenges are put forward by observations of 
magnetars. In Sect. [2] we briefly overview recent work on 
the construction of models of neutron star atmospheres 
with strong magnetic fields and on the modeling of spec- 
tra of thermal radiation formed in an atmosphere or at a 
condensed surface. We list important unsolved theoret- 
ical problems which arise in this modeling. In Sect. [Sj 
after a brief review of the effects of strong magnetic 
fields on the thermal structure and effective tempera- 
ture of neutron stars, we describe our new calculations 
of the thermal structure. Compa. r ed to the pre vious re- 
sults (jPotekhin fc Yakovlevll200l iPotekhin et al. 2003), 
we have taken into account neutrino energy losses in the 
outer crust of the star. We show that neutrino emission 
strongly affects the temperature profile in a sufficiently 
hot neutron star and places an upper limit on its surface 
temperature Tg and photon luminosity L^. 



2 Thermal Emission from Magnetized Surface 
Layers of a Neutron Star 



2.1 Neutron Star Atmosphere 



It was realized long ago ( Pavlov et al.lll995l ) that a neu- 
tron star atmosphere model should properly include the 
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effects of a strong magnetic field and partial ioniza- 
tion. Models of fully ionized neutron star atmospheres 
with strong magnetic fields were constructed by several 
research groups (e.g. , IShibanov et al.l Il992l : IZane et al] 
and references therein). The most 
recent papers highlighted the effects that can be impor- 
tant for atmospheres o f magnetars : the ion cyclotron 
feature (|Ho fc Lail[200ll: IZane et al .l l200lh and vacuum 
polarization, including a conversion of normal radia- 
tion modes propagating in the magnetized atmosphere 
( Ho k Lail[200alLai fc H0II2OO3D . 

Early studies of partial ionization in the magnetized 
neutro n star atmospheres (e.g., Rajagopal, Rom ani, & 
Miller Il997t reviewed by IZavlin &: Pavlovl |2002[ ) were 
based on an oversimplified treatment of atomic physics 
and nonideal plasma effects in strong magnetic fields. 
At typical parameters, the effects of thermal motion of 
bound species are important. So far these effects have 
been taken into account only for hydrogen plasmas. 
Thermodynamic functions, absorption coefficients, the 
dielectric tensor and polarization vectors of normal radi- 
ation modes in a strongly magnetized, partially ionized 
hydrogen plasma have been obtained and used to cal- 
cula te radiative opacities and thermal radiation spectra 
(see iPotekhin et al]l2004 and references therein). 

The summary of the magnetic hydrog en atmosphere 
model s and the list of references is given bv lPotekhin et all 
( 200d ). The model is sufficiently rehable at lO^^ G<B< 
2^013.5 j ^]-^g field range typical of isolated radio 
pulsars. It provides realis t ic sp ec tra of thermal X-ray 
radiat ion ( Potekhin et al ] l2004f ). iPotekhin fc Chabried 
(2003) extended this model to higher B. However, there 
remain the following unsolved theoretical problems that 
prevent to obtain reliable results beyond the indicated 
field range. 

— The calculated spectra a.t B > 10^^ G depend on the 
adopted model of mode conversion owing to the vac- 
uum resonance and on the description of the propa- 
gation of photons with frequencies below the plasma 
frequency. Neither of these problems has been def- 
initely solved. Their solution is also important for 
modeling the low-frequency (UV and optical) tail of 
the spectrum. 

— At low T or high B, hydrogen atoms recombine in 
H„ molecules and eventually form a condensed phase 
(see Sect. 12. 2p . Corresponding quantum-mechanical 
data are very incomplete. 

— At 10^G<i?< 10^^ G, transition rates of moving 
H atoms have not been calculated because of their 
complexity. There is the only one calculation of the 
energy spec trum of bound states ap propriate to this 
range of B (jLozovik fc Volkovll2004l ). 

— A more rigorous treatment of radiative transfer in the 
atmosphere requires solving the transfer equations for 
the Stokes parameters which has not been done so far 
for partially ionized atmospheres (see, e.g.. lLai fc Hoi 



I2OO3I : Ivan Adelsberg fc Lail[2006l for the cases of fully 
ionized atmospheres). 

Finally, we note that it is still not possible to cal- 
culate accurate atmospheric spectra at B > 10^^ G 
for chemical elements other than hydrogen, because of 
the importance of the effects of motion of atomic nu- 
clei in the strong magnetic fields. Apart from the H 
atom, the se effects have been calculate d only for the 
He atom (|A1-Huiai fc Schmelcheijl2003aU bh. which rests 
as a wh ole, but has a moving nucleus, and for the 
He+ ion ( Bezchastnov. Pavlov, fc Venturalll998l : Pavlov 
& Bezchastnov 2005( ). The data of astrophvsical rele- 
vance for He"*" are partly published and partly in prepa- 
ration (see iPavlov fc Bezchastnovl[2005l ): one expects to 
have a He/He"*" magnetic atmosphere model available in 
the near future. 



2.2 Condensed Surface and Thin Atmosphere 

The notion that an isolated magnetic neutr on star has 
a cond ensed surface was first put forwar d bvlRuderman 
1971 ). who considered the iron surface. ILai fc Salpeter 



1997t ) and 113 (|200ll ) studied the phase dia gram of 



strongly magnetized hydrogen and showed that, when 
the surface temperature Tg falls below some critical value 
(dependent of B), the atmosphere can undergo a phase 
transition into a condensed state. A similar phase tran- 
sition occurs for the equation of state of partially ion- 
ized, nonideal, strongly magnetized hydrogen p lasma , 
constructed by Potekhin, Chabri er, fc Saumon (Il999l) 
for B < 10^^-^ G and extended bv lPotekhin fc Chabrieil 
(I2OO4D ;o the magnetar field strengths. It is analogous 
to the "plasma phase transition" suggested in plasma 
physics a t B = (see, e.g., Chabrier, Saumon, fc 
Potekhin 2006 for discussion and references). Accord- 
ing to Potekhin et al. (j 19991) . the critical point for the 
phase transition in the hydrogen plasma is located at 
the density pc ~ 14:3Bl2^ g cm~'^ and temperature 
Tc « 3 X 10^B^2^^ K, where B12 = B/10^'^ G. At T < 
the density pcond of the condensed phase increases up 
to a few t imes of pc- On the other hand, according to 
Lail ( 200 ll ) , the surface density of a condensed phase for 
heavy elements is Pcond ~ 560AZ'^'^-^Bl2 g cm^"^, where 
Z and A are the charge and mass numbers of the ions. 
These two estimates of Pcond are in qualitative agree- 
ment. iLail ( 200ll ) estimates the critical temperature of 
hydrogen as < O.IE^ < lO^-^B^g* K {E^ being the 
cohesi ve energy), a lso in agreement with the above esti- 
mate. I.TonesI (|l98fih calculated the cohesive energy for 
Ne and Fe at lO^^ G< B < lO^^ G, using the den- 
sity functional theory (DFT ) and obtained Eg ~ 0.1i3i2 
keV. Recently, iMedin fc Lail (|2006h performed DFT cal- 
culations of the cohesive energies for zero-pressure con- 
densed hydrogen, helium, carbon, and iron at 10^^ G 
< B < 10^^ G. For instance, they found that the co- 
hesive energy per carbon atom ranges from ^ 50 eV at 
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B = 10^2 Q to 20 keV at lO^^ G. The cohesive energy 
per iron atom varies from ~ 0.8 keV at i? = 10^'^ G to 33 
keV at 10^^ G. These calculations suggest Tc of the same 
order of magnitude as the a bove estimate fo r hydro gen. 

Note that the mo dels of iPotekhin et al" and 
iLai fc Salpeted ( 1997t) are constructed in the framework 
of the "chemical picture" of plasmas, whose vahdity 
near the plasm a phase transition can be questionable 
(jChabrier et al.l [20001 . Thus the position (and the very 
existence) of the condensed surface requires further the- 
oretical investigation and experimental or observational 
verification. Hopefully, this can be done by analyzing ob- 
servations of thermal emission from neutron stars. 

The thermal emission from the magnetized sur- 
face was studie d by jB rinkm arml (Il980h. Turolla , Zane , 
& Drake (l2004l). iPerez-Azorm. Miralles. fc Pons! (|2005f ). 
and Ivan Adelsberg et al.l ( 2005 ). The spectrum exhibits 
modest deviations from blackbody across a wide en- 
ergy range, and shows mild absorption features associ- 
ated with the ion cyclotron frequency (energy fkUci = 
6.3B12Z/A eV) and the electron plasma frequency (en- 
ergy hujp — 2%^ pZjA eV, where p is in g cm^'^). How- 
ever, the predictions of the ion cyclotron feature and 
the spectrum at lower frequencies are not firm. The un- 
certainty arises from motion of the ions in the electro- 
magnetic field around their equilibrium lattice positions. 
Most of the models treat the ioii s as fixed (non- moving). 
Only lvan Adelsberg efaH (|2005f ) considered two limits of 
fixed and free ions. In reality, howe ver, the ions are nei- 
ther fixed nor completely free (see Ivan Adelsberg et al.l 
I2OO5I for estimates of possible uncertainties). 

In addition, the condensed surface of a neutron star 
can be surrounded by a "thin" atmosphere, which is 
transparent to X-rays, but optically thick at lower wave- 
lengths. Such a hypothesis ha s been first invoked by 
iMotch. Zavlin. fc Haberll ( 2003| l for explaining the spec- 
trum of the isolated neutron star RX J0720.4— 3125. 
Recently, the hydrogen atmosphere model, described in 
Sect. 12.11 together with the c o ndens ed surface emission 
model o f Ivan Adelsberg et all (|200i) have been success- 
fully used for fitting the spectrum of the isolated neu- 
tron star RX J1856. 5-3754 (H o et al . 2006) (assuming 
the atmosphere to be "thin" as defined above). 



3 Heat Transport through Magnetized 
Envelopes 

3.1 Overview of Previous Work 

The link between the magnetized atmosphere and stellar 
interior is provided by a heat blanketing (insulating) en- 
velope. The solution of heat transport problem relates the 
effective surface temperature Tg to the temperature Tb 
at the inner boundary of the blanketing envelope. With- 
out a magnetic field, it is conventional to place the inner 
boundary at p = pb 10^° g cm^'^. In this case, it can be 



treated in the quasi- Newtonian approximation with frac- 
tiona l errors < 10""^ ( Gudmundsson. Pethick. fc EpsteinI 
Il983f ). A strong magnetic field, however, greatly affects 
heat transport and, consequently, the thermal struc- 
ture of the envelope. The thermal structure of neutron 
star envelopes with radial magnetic f i elds ( normal to 
the surface) was studied bv IVan Riper] ( 19881 ) (also see 
IVan Riperj ll98^ for references to earlier work) . His prin- 
cipal conclusion was that the field reduces the thermal 
insulation of the heat blanketing envelope due to the 
Landau quantization of electron motion. The thermal 
structure of the envelope with magnetic fiel ds normal 
and ta nge ntial to the surf ace was analyzed bv iHernquisd 
(|1985D andlSchaaJ l|l990ah . The tangential field increases 
the thermal insulation of the envelope, because the Lar- 
mor rotation of electrons reduces the transverse electron 
thermal conductivity. 

The case of arbitrary angle 9b between the field 
lines and the normal to the surface was studied by 
iGreenstein fc Hartkd ( 19831 ) in the approximation of con- 
stant (density and temperature independent) longitudi- 
nal and transverse thermal conductivities. The authors 
proposed a simple formula which expresses Tg at arbi- 
trary 6b through two values of Tg calculated at = 
and 90 ° . The case of arbitrary 9b was stu died also by 
ISchaaj (fT990bl) andlHevl fc HernauistI ((199 8. 2001). 

iPotekhin fc YakovlevI (|200lh reconsidered the ther- 
mal structure of blanketing iron envel opes for any 9b, 
using improved thermal conductivities (jPotekh in 1999). 
IPotekhin et al.l ( 20031 ) analyzed accreted blanketing en- 
velopes composed of ligh t elements . In ag reement with 
an earlier conjecture of Hernauist (Il985) and simpli- 
fied treatments of lPagd (119951) and lShibanov fc YakovlevI 
([19 9 6) , they demonstrated that the dipole magnetic field 
(unlike the radial one) does not necessarily increase the 
total stellar luminosity at a given Tb. On the con- 
trary, the field B ~ 10^^-10^^ G lowers L y, and only the 
fields B > 10^^ G significantly increase it. IPotekhin et all 
(|2003D shifted the inner boundary of the blanketing enve- 
lope to the neutron drip, pb = 4 x 10^^ g cm""^, because 
in some cases they found a non-negligible temperature 
drop at p > 10^° g cm~^. They obtained that magnetized 
accreted envelopes are generally more heat-transparent 
than non-accreted ones (the sam e is true in the field-free 
case, studied by Potekhin, Cha brier. fc Yakovl cv 1997). 
However, this heat transparency enhancement is less sig- 
nificant, when the transparency is already enhanced by 
a superstrong magnetic field. 

Recentlv lPotekhin. Urpin. fc Chabrieil ( 20051 ) showed 



that qualitatively the same dependence of on B and 
on the chemical composition holds not only for dipole, 
but also for small-scale field configurations. 

iGeppert. Kiiker. fc Pagel (|2004L I2006D studied heat- 
blanketing envelopes with a magnetic field anchored 
either in the core or in the inner crust of the star 
(with dipole and toroidal field components of different 
strengths). They showed that a superstrong field in the 
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inner crust of a not too hot star can significantly af- 
fect the surface temperature distribution and make it 
nonuniform and even asymmetric, with hot spots hav- 
in g different temperatures. Similar result s were obtained 
by lPerez-Azorm. Miralles. fc Ponj ( 20061 ) , who also eval- 
uated pulsed fractions and phase-dependent spectra of 
neutron stars with strong magnetic fields anchored in 
the inner crust. 



3.2 Heat-Blanketing Envelopes of Magnetars: The 
Effect of Neutrino Emission 

Magnetars differ from ordinary pulsars in two respects: 
they possess superstrong surface magnetic fields, and 
they are generally younger and hotter. The first circum- 
stance suggests to extend the heat-blanketing layer to 
deeper layer s (at least to the ne utron drip density as 
was done bv'Potekh in et ani2003( ). The second indicates 
that neutrino emission can be important in the heat- 
blanketing envelopes. Accordingly, we have modified our 
computer code, w hich calculates the t herma l struc ture 
(|Potekhin k Yako vlev 2001; Potckhin'^aD 120031 ). to 
make it fully relativistic and to allow for energy sinks 
during heat diffusion. 



3.2.1 Basic Equations 

A complete set of equations for mechanical and thermal 
structure of a spherically sym metric star in h ydrostatic 
equilibrium has been derived bv lThorn^ ( 19771 ). They can 
easily be transformed to the form valid in an envelope 
of a star with radial heat transport, anisotropic (slowly 
varying) temperature distribution over any spherical 
layer, and a force-free magnetic field. Assuming quasis- 
tationary heat transport and neutrino emission, these 
equations reduce to the following system of ordinary dif- 
ferential equations for the metric function (gravitational 
potential) "P, the local heat flux Fr, temperature T, and 
gravitational mass m, contained within a sphere of cir- 
cumferential radius r, as functions of pressure P: 



dlnP 
1 d{r^Fr 



1 P 
P Q 



r2 dlnP 



dlnT 
dlnP 

dr 
dlnP 

dm 
dlnP 



Fr K'P 



- 2Fr 

1 



d<p 

dlnP' 
d<? 



16 crT^ g 

P ICr 

pg IChlCg ' 
Aur^PJCr 
glChJCg 



JChJCn dlnP' 



(1) 
(2) 
(3) 
(4) 
(5) 



Here p is the mass density (equivalent energy density of 
the matter), Q is the net energy loss per unit volume 
{Q = Qu is the neutrino emissivity in our case, although 



generally Q — — Qh, Qh being the heat deposition 
rate, e.g., due to nuclear reactions), K' = Kpba.i/p, K is 
the opacity, pbar = ^bar'TiH is the so-called "baryon mass 
density," ribar is the baryon number density, mpj is the 
mass of the hydrogen atom, a is the Stefan-Boltzmann 
constant, g = Gm/{r^lCr) is the local gravity, and G is 
the gravitational constant. Furthermore, 

/C, = (1 - 2Gm/rc2)i/2^ (6) 

ICh = l + P/pc\ (7) 



/Cg = 1 + inr^P/mc^ 



(8) 



are general relativistic corrections to radius, enthalpy 
and surface gravity, respectively. We adopt the conven- 
tional definition of the opacity K, used also by iThornd 
( 1977t ). In our notations K' = 16aT^/{3Kp), where k is 
the total (electron plus radiative) thermal conductivity 
of the plasma. Note that p « pbar and K' » K in the 
entire neutron star envelope. We use the same thermal 
conductivities as in iPotekhin fc YakovlevI (|200lh . The ef- 
fective radial thermal conductivity in a local part of the 
surface equals k = Ky cos^ 9b + k± sin^ 9b, where Ky and 
K± are the components of the conductivity tensor re- 
sponsible for heat transport along and across field lines, 
respectively. 

The local (non-neutrino) luminosity equals the inte- 
gral of the flux over the sphere of radius r. 



Lr 



sin 9d9 dipr^Fri9,ip), 



(9) 



where {9, (p) are the polar and azimuthal angles. Fo r 
a magnetic dipole model ( Ginzburg fc Ozernovl [196^ . 
tan0 — 2 tan^B. 

The boundary conditions to Eqs. H]) ~ Q are 

<Ps = In/Crs, Frs = Fr = aT^ , Ts = R, TOs = M, 

and the surface pressure is determined, following Gud- 
mundsson, Pethick, fc Epstein ( 19831 ). by the condition 
-ft'rad,s-Ps/5s = 2/3, whcrc iCrad is the radiative opacity. 
The subscript 's' refers to surface values; M is the grav- 
itational stellar mass. 

General Relativity correction factors ICr, K-h, and 
K,g in Eqs. (O - ([5]) are nearly constant because (A/ — 
m) /M ~ 10"^ and P/pc? ~ 10"^ at the bottom of the 
outer crust. However, they are taken into account in our 
code, in order to extend calculations to deeper neutron 
star layers, when required. 

Equations ([2]), ([3]) are one-dimensional, which implies 
that the mean temperature gradient along stellar radius 
is large compared to the tangential temperature gradi- 
ent, i.e., e = \dT /dx\/\dT / dr\ <^ 1, where a; is a coor- 
dinate along the stellar surface. Let us roughly estimate 
the mean value of e for a large-scale (e.g. , dipo le) mag- 
netic field, following I Greenstein fc Hartkd (|l983| l. In this 
case, e ^ {Ts/To){lo/R), where Iq is the depth at which 
the temperature distribution becomes nearly isotropic 
(T = To at the depth Iq) and Tg is the mean surface 
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logp [g cm-=] log/0 [g cm-3] 

Fig. 1 Temperature profiles in the outer crust of a neutron star with magnetic field B = 10^^ G (left panel) or _B = 10^® G 
(right panel), directed perpendicular to the stellar surface {9b = 0). Solid lines (NV) - ground-state (Negele-Vautherin) 
nuclear composition of the envelope, dot-dashed lines - ®^Fe envelope. For comparison, dashed lines show the temperature 
profiles without allowance for neutrino emission (for ground-state matter). For every family of curves, temperature at the 
neutron drip density is fixed to Tb = 10* K, 10^ K, and 2 x 10® K. The dots at the left end of the profiles correspond to the 
radiative surface, where the optical depth equals 2/3 and crT* = Fr. 



temperature. At the bottom of the outer crust (i.e., as- 
suming Tq — Tb) we have Iq/R ^ 0.1 (e.g., ~ 0-6 km 
for M = 1. 4A/0 and i? = 10 km) and Ts/Tq - 10"^ (ggg 
Sect. 13. 2. 2p . This estimate gives e < 10"^. 

Nevertheless, the one-dimensional approximation is 
inaccurate at those loci where magnetic field lines are 
tangential to the surface, because in this case one can- 
not neglect their curvature, dOs/dx. For large-scale mag- 
netic field {dOB/dx ~ the maximum size a of 
such sites can be estimated as the distance at which 
an initially tangential field line crosses the depth Iq 
(which assumes that heat flows along field lines, i.e., 
^ la/a). Thus, a < ^/ RIq. Since Tg is minimal 
at such sites, their contribution to the total stellar lumi- 
nosity can be neglected in the first approximation. 

3.2.2 Results 

We have solved Eqs. (H]) - ((51) using a straightforward 
generalization of the Runge-Kutta rnethod employed in 
our previ ous papers (iPotekhin et al. 1 119971 : Potekhin & 
Yakovlev 120011: IPotekhin et al.ll2003 [). Temperature pro- 
files have been calculated within a local part of the blan- 
keting envelope with a locally constant magnetic field. At 
every value of P, a corresponding value of p was found 
from the equation of s tate of magnetized electron -ion 
relativistic plasma (e.g. . IPotekhin k, Yakovlevll 19961 ). us- 
in g approximations of Fermi-D irac integrals presented 
in IPotekhin fc Chabrierl ( 2000f ). In all examples shown 
in the figures we chose the "canonical" neutron star 
mass M = I.AMq and radius i? = 10 km. The neu- 
trino emissivity is calculated as Qi, — Qpair + Qpi + 



Qsyn + Qbrcms, whcrc the Contributions due to electron- 
positron pair annihilation Qpair, plasmon decay Qpi, syn- 
chrotron radiation of neutrino pairs by electrons Qsyn, 
and bremsstrahlung in electron-nucleus collisions Qbrems 
ar e given, respecti v ely, by Eqs. (22), (38), (56), and (76) 
of Yakovlev et al. ('200 lJ). According to th e results of 
lEoh et al. (1996) and Yakovlev et all (|200ll ). other neu- 
trino emission mechanisms are unimportant in the outer 
crust of neutron stars. 

Our calculations show that neutrino emission is cru- 
cially important for the thermal structure of neutron 
stars with internal temperature Tb > 10^ K. 

Figure [T] shows temperature profiles at Tb = 10* K, 
10^ K, and 2 X 10^ K for magnetic fields B = 10^^ q 
and 10^^ G, directed perpendicular to the stellar sur- 
face. The present results (solid and dot-dashed lines) are 
compared with the profiles calculated neglecting neu- 
trino emission (dashed lines). At the lowest tempera- 
ture Tb — 10* K there is virtually no difference (all the 
lines coincide). At Tb — 10^ K, the difference is notice- 
able, and at Tb = 2 X 10® K it is large. Because of the 
growth of neutrino emission with increasing temperature 
at p > 10^" g cm~'^, Tg is nearly independent of Tb at 
Th }t 10® K, but depends on the magnetic field. 

In this figure we also compare temperature profiles for 
the different heavy-element compositions of the outer en- 
yelope: iron (dot-dashed li nes) and ground-state matter 
( Negele fc VautheriIi^ll973^ . The effect of composition is 
not strong, but noticeable when the neutrino emission is 
important, because Qi, depends on the electron number 
density that is a function of composition for a given p. 
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Fig. 2 Profiles of local radial heat flux Fr for the cases shown 
in Fig.[l] Top panel: B = 10^^ G, bottom panel: B = 10^^ G; 
solid lines - ground-state matter, dot-dashed lines - ^^Fe, 
dashed lines - without neutrino emission for ground-state 
matter. 



Figure demonstrates the profile of the local heat 
flux Fr, for the same cases as in Fig. [1] plotted by 
the same line types. Without neutrino emission (dashed 
lines), Fr would be nearly constant, with only w 2% in- 
crease towards the inner crust due to the General Rela- 
tivity effects (associated with the variation of the metric 
function <P) and « 9% increase because of the spherical 
geometry (the factor). The neutrino emission leads to 
a strong dependence of the flux on p and violates the fa- 
miliar relation between Fr and Tg derived in the absence 
of energy sinks. 

Figure [3] shows which neutrino emission mechanism 
dominates at given p, Tb, and B. For superstrong mag- 
netic fields, the neutrino synchrotron mechanism domi- 
nates in certain density ranges, which does not happen 
at "ordinary pulsar" B ^ 10^^ G. It is natural because 
of the strong i?-dependence of the synchrotron neutrino 
emissivity. Pair annihilation neutrino emissivity is not 
seen in the figure, because it is too small. Notice that 
the emissivity of plasmon decay and bremsstrahlung pro- 
cesses can be affected by superstrong magnetic fields 




10 11 
logp [g cm-3] 



10 11 
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Fig. 3 Contributions to the total neutrino emissivity 
Qi, from plasmon decay (solid lines), electron neutrino 
bremsstrahlung (dot-dashed lines), and electron synchrotron 
radiation (dashed lines). Left panels: B = W^^ G (radial 
field); right panels: B = 10^^ G; top panels: T = 10^ K; 
bottom panels: T = 2 x 10'' K. 



which has not been explored so far (and we present the 
emissivities in the field-free case). A slow dependence of 
these emissivities on B, seen in Fig.[3l is indirect (caused 
by the dependence of temperature profiles on B). 

Figures m and O give the Ts(rb) relation for the mag- 
netic fields B — 10^^ G and B = 10^^ G perpendicular 
and parallel to the radial direction. The relations ob- 
tained with and without allowance for neutrino emission 
are plotted by solid and dashed lines, respectively. We 
see that at Tb < 10* K the neutrino emission does not 
affect Tg. At higher Tb > 10® K, in contrast, this emis- 
sion is crucial: if = 0, then Tg continues to grow up 
with increasing Tb, whereas with realistic the surface 
temperature tends to a constant limit, which depends 
on 9b and B. In most cases this limit is reached when 
Tb ^ 10® K, but for a superstrong field (right panel) and 
transverse heat propagation, it is reached at still smaller 
Tb ~ 3 X 10* K. 

In Fig. [5] we explore joint effects of magnetic field, 
neutrino emission, and the shift of the inner boundary 
from 10^'^ g cm~'^ to the neutron drip. Here magnetic 
field lines are directed along the surface (perpendicu- 
lar to the direction of heat transport). Therefore, these 
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Fig. 4 Surface temperature Ts as a function of tempera- 
ture Tb at the neutron drip point for a neutron star with 
the magnetic field B = 10^^ G, directed along stellar radius 
{6b = 0, parallel heat transport, sign ||) or along the surface 
{9b = 90°, transverse transport, _L). Ground-state composi- 
tion is assumed. Solid lines - present calculation, dashed lines 
- neutrino emission is neglected. 



effects are most pronounced. The solid lines are the re- 
sults of accurate calculations; the dashed lines, as before, 
show the model with — 0; and dot-dashed lines are 
obtained by solving Eqs. ([1]) - ([5]) in the whole domain 
Ps < P < 4 X 10^^ g cm"'^, but with the magnetic field 



artificially "switched off" at p > 10"'^'' g cm~'^. This is a 
simulation of the model, where the heat transport in the 
magnetized plasma is solved accurately up to pb = 10^" 
g cm^'^, while after this boundary nonmagnetic heat bal- 
ance and transport equations are solved. In the absence 
of the result s repor t ed he re, the latter model was used 
by [Kamin ker et al.l ( 2006t) to study the thermal struc- 
ture and evolution of magnetars as cooling neutron stars 
with a phenomenological heat source in a spherical in- 
ternal layer. In the left panel of Fig. [6l the solid and 
dot-dashed lines almost coincide, indicating that in this 
case ph = 10^" g cm~'^ may provide a sufficient accuracy. 
In the right panel, in contrast, the effect of the shift of the 
inner boundary is quite visible. Therefore, we conclude 
that the development of our thermal structure code, re- 
ported here, will allow us to study the thermal history 
of magnetars at a higher accuracy level. 

Meanwhile, a comparison of the profiles shown in 
Figs.[T]and[6]prompts that at T > 10^ K and J5 ~ lO^^ G 
the magnetic effects on the conductivity could be impor- 
tant at still higher densities in the inner crust. Investi- 
gation of this possibility requires taking into account the 
effects of free (possibly superfluid) neutrons on thermal 
conduction and neutrino emission. We are planning to 
perform such study in the future. 

In the above figures we have shown the results of 
calculations where the surface of a neutron star was 




8 8.5 
log (K) 

Fig. 5 The same as in Fig. g] but for B = 10^^ G. 



assumed diffuse, i.e., without the phase transition dis- 
cussed in Sect. [521 However, for B = 10^^ G and Tb < 
10^ K the surface temperature Tg is below the estimates 
of the critical temperature Tc given in Sect. 12.21 This pos- 
sibility is explored in Fig. [71 Here we use the equation of 
state containing nonideal terms fo r strongly magnetize d 
fully ionized plasma (Sect. IIIB of jPotekhin et al]|l999l ). 
which enforce the phase transition. At B = 10^^ G the 
surface density is very high, pcond ~ 3.1 x 10^ g cm~^, in 
our model (and the analytic estimate in Sect. 12.21 gives 
Pcond ~ 2 X 10^ g cm~^). The solid lines in the left panels 
reproduce the profiles shown in Figs. [J and [HI whereas 
the dot-dashed lines display the case of magnetic con- 
densation. In the right panels, we show the same tem- 
perature profiles as a function of local proper depth I 
{dl = — /C~^dr), measured from the radiative surface. Al- 
though the temperature profiles with and without mag- 
netic condensation are drastically different in the sur- 
face layers, the effective temperature remains almost the 
same. Thus, the magnetic condensation does not signif- 
icantly affect the Ts(Tb) relation. This should not be 
surpri sing, because, as explained by lGudmundsson et al.l 
(I1983D . the main regulator of the Tb - Tg relation is the 
"sensitivity strip" where k has a minimum. This domain, 
a real bottleneck for the heat outflow, lies at p > Pcond 
(except for low Tg, not considered here), and therefore it 
is not affected by the condensation. 



4 Summary 

We have described the main effects of strong magnetic 
fields on the properties of neutron star atmospheres and 
heat blanketing envelopes, with the emphasis on the dif- 
ference between ordinary neutron stars and magnetars. 
Observations of magnetars pose new theoretical prob- 
lems and challenges because of magnetars' supcrstrong 
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Fig. 6 Temperature profiles in tlie outer crust of a neutron star witli tlie magnetic field B — 10^^ G (left panel) or B = 10^"' G 
(right panel), directed parallel to the stellar surface, with the temperature Tb = 10^ K at neutron drip density and the ground- 
state composition of the matter. Solid lines - present calculation, dot-dashed lines - with a switch to non-magnetic {B — 0) 
calculation at p > 10^'^ g cm~^, dashed lines - without neutrino energy losses. 



magnetic fields and high temperatures. We also report 
a solution to one of such problems, which consists in 
taking into account neutrino energy losses in the outer 
crust of hot, strongly magnetized neutron stars. We have 
demonstrated that, because of these losses, the effective 
surface temperature Tg almost ceases to depend on the 
temperature Tb in the inner crust as soon as Tb exceeds 
10^ K. A direct consequence of this observation is that 
in the absence of powerful energy sources in outer en- 
velopes, the stationary (time-averaged) effective temper- 
ature cannot be raised above the value that it would have 
at Tb « 10^ K, irrespectively of the energy release in the 
deeper layers. 
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